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Abstract 

The effective theory for baryons with combined l/N^ and chiral expansions is analyzed for non- 
strange baryons. Results for baryon masses and axial couplings are obtained in the small scale 
expansion, to be coined as the ^-expansion, in which the and the low energy power countings 
are linked according to 1/Nc = 0{Cj = 0{p). Masses and axial couplings are analyzed to C'(^^) 
and 0{^^) respectively, which correspond to next-to-next to leading order evaluations, and require 
one-loop contributions in the effective theory. The role of the spin-flavor approximate symmetry 
in baryons, consequence of the large Nc limit, is manifested in the physical world with Nc = 'i 
in a significant way, as the analysis of its breaking in the masses and the axial couplings shows. 
Applications to the recent lattice QCD results on baryon masses and the nucleon's axial coupling 
are presented. It is shown that those results are naturally described within the effective theory at 
the order considered in the (^-expansion. 
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I. INTRODUCTION 

The low energy effective theory for baryons is a topic that has evolved over time through 
several approaches and improvements. The early version of baryon Chiral Perturbation 
Theory (ChPT) [1] evolved into the various effective field theories based on effective chiral 
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Lagrangians [2-4], starting with the relativistic version [5, 6] or Baryon ChPT (BChPT), 
followed by the non-relativistic version based in an expansion in the inverse baryon mass 
[6, 7] or Heavy Baryon ChPT (HBChPT), and by manifestly Lorentz covariant versions based 
on the IR regularization scheme [8-10]. In all these versions of the baryon effective theory 
a consistent low energy expansion can be implemented. The most important issue, which 
became apparent quite early, was the convergence of the low energy expansion. Being an 
expansion that progresses in steps of 0{p) in contrast to the expansion in the pure Goldstone 
Boson sector where the steps are O(p^), it is natural to expect a slower rate of convergence. 
However, a key factor with the convergence has to do with the important effects due to 
the closeness in mass of the spin 3/2 baryons. It was realized [11], that the inclusion of 
those degrees of freedom play an important role in improving the convergence of the one- 
loop contributions to certain observables such as the vr-A^ scattering amplitude and the axial 
currents and magnetic moments. There have been since then numerous works including spin 
3/2 baryons [12-20]. The key enlightenment resulted from the study of baryons in the large 
Nc limit of QCD [21]. It was realized that in that limit baryons behave very differently than 
mesons [22], in particular because their masses scale like 0{Nc) and the vr-baryon couplings 
are 0{\/Nc). These properties were shown to require for consistency that at large Nc baryons 
must respect a dynamical contracted spin flavor symmetry SU{2Nf), Nf being the number 
of light flavors [23-26], broken by effects ordered in powers of l/N^ and in the quark mass 
differences. The inclusion of the consistency requirements of the large Nc limit into the 
effective theory came naturally through a combination of the 1/Nc expansion and HBChPT 
[27], which is the framework followed in the present work. The study of one-loop corrections 
in that framework was first carried out in Refs. [27-29]. In the combined theory one has to 
deal with the fact that the 1/Nc and Chiral expansions do not commute [30]. The reason is 
due to the presence of the baryon mass splitting scale of 0{1/Nc) {A — N mass difference), 
for which it becomes necessary to specify its order in the low energy expansion. Thus the 
1/Nc and Chiral expansions must be linked. Particular emphasis will be given to the specific 
linking in which the baryon mass splitting is taken to be 0{p) in the Chiral expansion, and 
which will be called the ^-expansion. Following references [27-29], the theoretical framework 
is presented here in detail, in particular the renormalization, the power countings, and the 
linked 1/Nc and low energy expansions, along with observations that further clarify the 
significance of the framework. 
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The very significant contemporary progress in the calculations of baryon observables in 
lattice QCD (LQCD) [31-33] opens new opportunities for further understanding the low 
energy effective theory of baryons. The determination of the quark mass dependence of the 
various low energy observables, such as masses, axial couplings, magnetic moments, elec- 
tromagnetic polarizabilities, etc., are of key importance as a significant test of the effective 
theory, in particular its range of validity in quark masses, as well as for the determination 
of its low energy constants (LECs). Lattice results for the and A masses [34-39] and 
the axial coupling of the nucleon [40-45] at varying quark masses are analyzed with the 
purpose of testing the effective theory presented here. This in turn can give insights on 
LQCD results, in particular an understanding on the role and relevance of including the 
spin 3/2 baryons consistently with large Nc requirements. 

This work is organized as follows. In Section II the framework for the combined 1/Nc 
and HBChPT expansions is presented. Section III presents the evaluation of the baryon 
masses and Section IV the one for axial couplings at the one-loop level. Section V is de- 
voted to applying those results in the ^-expansion to LQCD results. Finally, Section VI is 
devoted to observations and conclusions . Several appendices present useful material used 
in the calculations, namely. Appendix A on spin-flavor algebra. Appendix B on symmetries, 
Appendix C on the construction of effective Lagrangians, and Appendix D on useful matrix 
elements of spin-flavor operators. 



II. FRAMEWORK FOR THE COMBINED 1/Nc EXPANSION AND BARYON 
CHIRAL PERTURBATION THEORY 

In this section the framework for the combined 1/Nc and chiral expansions in baryons is 
presented in some detail along similar lines as in the original works [27-29]. The symmetries 
that the effective theory must respect in the chiral and large Nc limits are chiral SUiiNf) x 
SUii{Nf) and contracted dynamical spin-flavor symmetry SU{2Nf)[23-2Q] ^. Nf is the 
number of light flavors, and in this work Nf = 2. In the limit A^^c — t- oo the spin-flavor 
symmetry requires baryons to belong into degenerate multiplets of SU{A). In particular, the 
ground state (GS) baryons belong into a symmetric SU{A) multiplet, which consists of states 

^ See also Appendix B. 
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with I = S, being 5* the baryon spin and / its isospin. At finite Nc the spin-fiavor symmetry 
is broken by effects suppressed by powers of l/Nc, and the baryon mass sphttings in the GS 
multiplet are proportional to {S + 1)/Nc. The effects of finite are then implemented as 
an expansion in 1/Nc at the level of the effective Lagrangian. Because baryon masses scale 
as proportional to Nc, it becomes natural to use the framework of HBChPT [7, 46], where 
the expansion in inverse powers of the baryon mass becomes part of the 1/Nc expansion. 
The framework presented next follows that of Refs. [27, 28]. 

The non-relativistic baryon field, denoted by B, consists of the symmetric spin-fiavor 
SU{4:) multiplet with states I = S, S = 1/2, ■■■ ,Nc/2 {N^ odd). Chiral symmetry is 
realized in the usual non-linear way on B, namely [2-4]: 

{L,R) -.B = h{L,R,u)B, (1) 

where L(R) is a SUl{r){'2) transformation, u is given in terms of the pion fields tt" by 
u = exp{iTT"'T°' /2Ft^), = 92.4 MeV, and h{L,R,u) is an SUi{2) isospin transformation 
which in any representation of Isospin satisfies Ruh'^{L, R, u) = h{L, R, u) u L'^ . The chiral 
covariant derivative -D^B is given by: 

D,B = d^B - iTltB, 

= ^(^^r,), (2) 

where = v^ — and = v^ + a^ are gauge sources. Here the notation (A) = TiA is used 
for fiavor traces, and in general. A"' = ^{t°'A) where A is in the fundamental representation 
in the trace. The definition = I2x2 is used. The axial Maurer-Cartan one-form necessary 
as building block of the effective chiral Lagrangian is: 

u^ = u\id^ + r^)u-u{id^ + l^)u\ {L, R) : u^, = h{L, R,u)Uf,h\L, R,u). (3) 

Since = 0{\^Nc), u, u^, contain different orders in the expansion in powers of 1/Nc. 
The contracted S'f/(4) transformations (see Appendix A) are generated by {5**, X*"}, 
where = G^'^ /Nc are semiclassical at large A^^- The ordering in Nc of the matrix elements 
of the spin-flavor generators in states with S = 0{N^) are as follows: S** = 0{N^), = 
0{N^), and G*" = 0{Nc). While inflnitesimal SU{4) transformations generated by J" 
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correspond to the usual isospin transformations when acting on pions, the ones generated 
by X*" affect only the baryons (one can define these generators to not affect the pion field as 
shown in Appendix B). The effective Lagrangian can be systematically written as a power 
series in the low energy expansion or Chiral expansion, and simultaneously in l/N^. It is 
most convenient to write the Lagrangian to be manifestly chiral invariant as is usually done. 
The low energy constants (LECs) will themselves admit an expansion in powers of 1/Nc- 
For the HBChPT expansion the large mass of the expansion is taken to be the spin-flavor 
singlet component of the baryon masses, Mq = NcmQ (mo can be considered here to be a 
LEG defined in the chiral limit and which will have itself an expansion in l/N^). To 0{1/Nc) 
baryon masses will read [25, 26]: 

m^{S)=M, + ^S{S + l) + ciN,Ml + --- . (4) 

The baryon mass splittings due to the hyperfine term, second term in Eq. (4), must be 
considered to be a small energy scale. It becomes necessary to establish of what order that 
term is in the low energy expansion, as it naturally appears in combinations with powers 
of when loop diagrams are calculated. This fact makes that the low energy and 1/Nc 
expansions do not commute [30, 47], and the natural way to proceed is therefore to link the 
two expansions. For the purpose of organizing the effective Lagrangian it is convenient to 
established the link between the two expansions. In the real world with A^^c = 3 the A — N 
mass splitting is about 300 MeV, and therefore it is reasonable to count that quantity as 
0{p) in the low energy expansion: the expansion where = C>{p) = 0{^) will be adopted 
in what follows, and it will be called ^-expansion. This power counting corresponds to the so 
called small scale expansion (SSE) [13], now consistently implemented in the context of the 
1/Nc expansion. Whenever appropriate, it will be indicated which aspects of the analysis 
are general and which are only valid in that expansion. Up to 0{C,) the baryon effective 
Lagrangian reads: 

= Bt (^^Do + QAU^aG^'' - - |Ar, x^j B, (5) 

where qa is the axial coupling in the chiral and large Nc limits (it has to be rescaled by a 
factor 5/6 to coincide with the usual axial coupling as defined for the nucleon), x+ is the 
source containing the quark masses: specifically x+ = 2M^ + ■ ■ ■ (see Appendix C ). Here 
one notes an important point which will be present in other instances as well: the baryon 
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mass dependence on the current quark mass behaves at 0{Nc M^) (ci is of zeroth order in 
Nc)i and this indicates that in a strict large Nc hmit the expansion in the quark masses of 
certain quantities such as the baryon masses cannot be defined due to divergent coefficients of 
0{Nc). The Lagrangian is manifestly invariant under the chiral transformations, translations 
and rotations (the latter also involving obviously the action of the 5** generators of SU{A)). 
Under transformations generated by X*", the kinetic term changes at 0{1/N^), the term 
proportional to gA which contains the ttBB' interaction and the leading order terms of the 
axial currents, changes by terms which are a factor 0{1/N^) smaller than the original term, 
the term proportional to Ci, which gives the leading order (LO) a-term in the baryon masses, 
is a spin-flavor singlet, and flnally the hyperflne term proportional to Chf is the only one 
that shows explicit spin-flavor symmetry breaking as it changes by 0{1/Nc), which is the 
order of the term itself (this is so because [S"^, X*"] — 0{N^)). That term therefore provides 
the dominant spin-flavor symmetry breaking in the effective Lagrangian. The higher order 
Lagrangians can be built using the tools provided in Appendix C . The operators appearing 
in the effective Lagrangian are normalized in such a way that all the LECs are of zeroth 
order in Nc- The 1/Nc power vi/n^ of the operators in the Lagrangian are determined by 
the following simple rule: 

i^i/jVe ^^--l-z^+y, (6) 

where the spin-flavor operator is n-body {n is the number of factors of SU{^) generators 
appearing in the operator), k is basically the number of factors of the generators G*" remain- 
ing after reducing the operator using commutators, and n^^ is the number of pions attached 
to the vertex. It is opportune to point out that commutators of spin-flavor operators will 
always reduce the n-bodyness of the product of operators: e.g., let Q be any generator of 
SU{A), and consider the commutator [^,5"^] = {S'^ ,[Q , S'^]} , and then, because [Q,S'^\ is a 
1-body operator, [^,-5^] is actually a 2-body operator. 

A. Consistency of the 1/iVc expansion 

The consistency of the 1/Nc expansion in QCD gives rise to the dynamical spin-flavor 
contracted SU {2Nf) symmetry in baryons at large N(.. At the baryon level that symmetry 
can be deduced as the result of consistency or correct Nc power counting, of observables 
in which pion-baryon couplings are involved: because the pion-baryon coupling is 0{^/Nc) 
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from Witten's counting rules [22]. In particular the consistency of pion-baryon scattering 
is a direct way of deriving the existence of the dynamical spin-flavor symmetry [25, 26]. In 
general, for any quantity there must be cancellations between the terms with the "wrong" 
power counting stemming from different Feynman diagrams. For instance, baryon masses 
are 0{Nc), and therefore pion loop contributions cannot give contributions which scale with 
a higher power of Nc- On the other hand, the baryon mass splittings are 0{1/Nc), and 
loop contributions must respect that scaling. Similarly, in the axial currents, whose matrix 
elements are 0{Nc) such cancellations occur when loop corrections are calculated. All this 
will be illustrated in the application to baryon masses and axial couplings discussed next. 
Although certain key cancellations must be exact in the large Nc limit, the analysis of LQCD 
results will show that they are very significant in the physical world where Nc = 3. 

B. ^ power counting 

The terms in the effective Lagrangian are constrained in their Nc dependence by the 
requirement of the consistency of QCD at large Nc- This constraint is in the form of a lower 
bound in the power in 1/Nc for each term one could write down in the Lagrangian. This 
leads to constraints on the A^^^ dependencies of the ultra-violet (UV) divergencies, which have 
to be subtracted by the corresponding counter-terms in the Lagrangian. One very important 
point to mention is that the UV divergencies are necessarily polynomials in low momenta 
p (derivatives), in and in 1/A"c (modulo factors of I/a/Ac due to I/Ft^ factors in terms 
where pions are attached). Therefore, the structure of counter-terms is independent of any 
linking between the 1/Ac and chiral expansions. For this reason, one can simply take the 
large Nc and low energy limits independently in order to determine the UV divergencies. 
The Chiral or low energy order of a diagram is given by [48]: 



where ub is the total number of external baryon legs (two in the present work), L the number 
of loops, rii the number of vertices of type i, Up. is the the low energy (Chiral) power of that 
type of vertex, and bi the number of baryon legs of the vertex (0 or 2 in the single baryon 
sector). On the other hand, the 1/A"c power of a diagram is given by: 




(7) 
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where n.,^ is the number of external pions of vertex i, and the 1/Nc order of the spin- 
flavor operator of that vertex. Since z/o^ can be negative (due to factors of G*" in vertices), 
one can think of individual diagrams with Vi/Na negative and violating large N^. consistency, 
requiring cancellation with other diagrams. Such a sum will have to respect the mentioned 
lower bound on the 1/Nc power corresponding to the sum of such diagrams. The explicit 
example of such cancellation in the axial currents at one-loop is given below. One interesting 
point is that the leading in l/N^. order of terms in the Lagrangian are primarily determined 
by the order of the spin-flavor operators according to Eq.(6). 

One can determine now the nominal counting of the one-loop contributions to the baryon 
masses and axial currents. The one loop correction shown in Fig. 1 has : {L = l^ns = 
2,n7r = 0,ni = 2, i/Qi = — = 2,z/pj = 1) giving z/p = 3 as it is well known, and 
^i/Nc = Since there is only one possible diagram this must be consistent, and it is. For 
the axial currents one has the diagrams in Fig. IV. The current at tree level is (9(iVc), and 
the sum of the diagrams cannot scale like a higher power of A^^^. Performing the counting 
for the individual diagrams one obtains: t'p(j) = 2 for j = 1, ■ ■ ■ , 4, and ui/nXj) = ~2, 
j = 1,2,3 and z/i/7v^(4) = 0. Thus a cancellation must occur of the 0{N^) terms when the 
contributions to the axial currents by diagrams 1, 2 and 3 are added. Since the acceptable 
bound is that the sum be 0{Nc), one concludes that the axial current has, at one-loop, 
corrections 0{p'^N^ or higher. 

One can consider the case of two-loop diagrams. Consider diagrams where the same 
pion-baryon vertex Eq.(5) appears four times. For the masses one has t'p(j) = 5, and 
individual diagrams give vi/m^ = —2. A cancellation must occur to restore the bound on 
the Nc counting for the masses, i.e., 0{Nc). Thus, at two-loops the UV divergencies of the 
masses must be 0{p^Nc) or higher. For the axial currents a similar discussion requires that 
counter-terms to the axial currents must be 0{p'^N^ or higher. One loop with the term 
proportional to Ci in Eq.(5) give 0{^*) contributions to the masses. 

In order to determine the power counting in for a given Feynman diagram with a single 
baryon flowing through the diagram, it is sufficient to know the l/N^ power counting of the 
operators in the vertices, the number of external pions, and the chiral power counting of the 
vertices. Eliminating the numbers of pion and baryon propagators in the diagram using the 
known topological formulas [48]. A quick calculation then permits to express the order in ^ 
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of the diagram according to: 

i/^ = 1 + 3L + + ^ Uii^i +pi-l) (9) 

i 

where L is the number of loops, is the number of external pions, rij the number of vertices 
of type i, and t/j is the l/iV^ order of the spin-flavor operator appearing in and Pi the low 
energy (chiral) power of that type of vertex. It should be emphasized that because Uj, can 
be negative (due to Nc enhancement when the generator G'" is present as a factor in a 
vertex). The ^-power counting of the UV divergencies is obvious from the earlier discussion. 
At one-loop one finds that the masses have and 0{^^) counter-terms, while the axial 

currents will have 0{^) and 0{^^) counter-terms. To two loops one expects and 
0{^^), and 0{^^) and counter-terms for masses and axial currents respectively. The 

non-commutativity of limits is manifested in the finite terms where and or momenta 
and 5m appear combined in non-analytic terms, and are therefore sensitive to the linking of 
the two expansions. 

III. BARYON MASSES 

In this section baryon masses are analyzed to order or next-to-next to leading order 
(NNLO), in the limit of exact isospin symmetry. To that order the mass of the baryon of 
spin S reads: 

msiS) = N,mo + ^S{S + 1) + c^N^M^ + <5m^-''"'^+^^(5), (10) 

where Sm^''°°^'^^^ (S) involves contributions from the one-loop diagram in Fig. 1, and 
CT denotes counter-terms. Prom both types of contributions, there are C(^^) and 0(^^) 

terms, and the calculation is exact at the latter order, as can be deduced from the 
previous discussion on power counting. The notation for the 0{^) mass shift is used: 
Sm{S) = ^S{S + 1) + CiN^Ml. Notice that Chf is equal to the LO term in Ma - Mjv in 
the real world Nc — 3. 



The leading 1-loop correction to the baryon self energy, diagram in Fig. 1, and reads: 



S^i-ioop) = ^ || ^ ^G^'^KG- /(i-ioop)(5m„ -p°,M,) , (11) 
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FIG. 1: One- loop contribution to baryon self energy. The red propagator indicates sum over all 
possible baryons that can contribute. 

where n indicates the possible intermediate baryon spin-isospin states in the loop, Vn are 
the corresponding spin-flavor projection operators, 5mn = 6m{Sn), and the loop integral is 
calculated in dimensional regularization with the result, 
,^ , , d'^k P 1 

^{l~loop)[Q, Mt,) 



{2tiY k^-M^ + ie k^-Q + ie 

Q ml - 2Q')(Ae - log ^) + ml - 



+ 2^{Ml - Q^fl^ + 4(Q^ - ^')'^'tanh-^ -^=^=1 , (12) 

where Q = drrin — p^, \e = ^ — 7 + log47r, and /i is the renormalization scale which will 
be taken later to be of the order of nip. For the specific evaluation of 5S(i_ioop) for a given 
baryon state denoted by in, = 6min — j3°, where p° is a residual energy (when evaluated 
on an on-shell baryon it is the kinetic energy which is 0{p'^ /Nc)). The non-commutativity 
of the 1/Nc expansion and the Chiral expansion of course resides in the non-analytic terms 
of the loop integral though their dependence on the ratio Q^/M^. 

Appendix D provides all the necessary elements for the evaluation of the spin-flavor 
matrix elements in Eq. (11) as well as in the calculation of the one-loop corrections to the 
axial currents below. 

The one-loop contribution to the wave function renormalization constant is given by: 

d_ 

The explicit evaluation of the ultraviolet divergent pieces of the self energy gives: 

= (14) 



(13) 



X 1^ (^-3M2(3iV,(4 + AT,) - 205^) + 8^(iV,(4 + iV,)(3 + 5^^) - 45^(6 + 7S')) 

(jyf2 o (^2 \ ■) 

_E(^iV,(4 + iV,) - S') - ^(iVe(4 + iV,)(3 + 2S') - 8S'i3 + S'))j + 0{p'')j 
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The UV divergent pieces start at Note that the UV divergencies in the mass (term 

independent of p°) is produced by the contribution of the partner baryon and is proportional 
to the mass sphtting. As is well known they are absent in HBChPT without exphcit A. 
The 0{N^) UV divergence is spin-fiavor singlet and proportional to M^, while the contribu- 
tions to mass splittings are 0{1/N^). Notice that the leading UV divergence of 6Zi_ioop is 
0{M^Nc): this is necessary as shown later for rendering the one-loop calculation of the axial 
currents consistent in the large Nc limit. Since the calculation is accurate to addi- 
tional terms in the effective Lagrangian up to that order are necessary for renormalization. 
The terms necessary for renormalizing the self energy are therefore the following: 

where one identifies the residual energy p° with the operator {iDo — 5m). All LECs are here 
of the form X{Nc) = Xo + Xi/Nc + ■ ■ ■ ■ Writing X = X{fi) +7^ Ae, one rcnormalizes the self 
energy to 0{^^). The coefficients 7^ are determined from ^S^^^^^^^^ given above. While the 
counter-terms are defined such that X{ijl) is 0{N^), it is possible that 7^ is of higher order 
in l/Nc- Notice that among the higher order terms there are terms which can be simply 
absorbed into l/Nc corrections to the LECs of the lowest order Lagrangian, and to ttiq. 
Finally, the baryon masses are given by: 

mB = Kmo + + cKMl + {STF^I/J:;^^ + 5E^^) |,o=o (1 + ^Z^JiJ^^^'^ + 6Z^^) (16) 

Note that the correction to the wave function renormalization factor enters in the expression 
for the mass corrections: this is because (^S(p° — 0) starts with terms 0[^^) and 5Z starts 
at 0(0) therefore the 0{^^) terms of the mass correction involve these lower order terms of 
the wave function renormalization. 

The one-loop corrections and corresponding counter-terms contribute to the masses at 
0{^'^) and O(C^), while in a strict large Nc limit the following ordering is found: 

Mb = 0{Nc) + 0{NcMl) + 0{N^,M^) + • • • 
M.-M^ = C.(i-) + 0(^) (17) 

Obviously the term 0{^^^^) stems from the 1/A^c expansion of non-analytic terms and 
shows the non-commutativity of limits. 
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The (T-terms for and A, defined hy as = ^'^^^ = \ijrLu + ^^d)), with the 

one-loop corrections contributing up to 0{S,^). The difference — a^' is and at that 

order it receives only finite contributions from the loop. This implies that the slopes of the 

and A masses as functions of M^^ are the same up to deviations. This seems to be 

closely followed by the lattice QCD results analyzed later. In the large Nf. limit, obviously 
cr = 0{N^. The terms of that order are necessarily spin-flavor singlet, and taking the limit 
at fixed M^r one finds cta — <^n = 0{^/^c)-> ^ result similar to the one in the ^-expansion. 



IV. AXIAL COUPLINGS 

In this section the evaluation of the axial couplings including corrections O(C^) is pre- 
sented. At that order the one-loop corrections must be calculated. 

The matrix elements of interest for the axial currents are (B' | A''^ \ B) evaluated at 
vanishing external 3-momentum. The axial couplings are then defined by: 

(B' I A'"" I B) = g^^' - (B' I G''' \ B) . (18) 

6 

The axial couplings defined here are 0{N^). The 0{N^ of the matrix elements of the axial 
currents is due to the operator G*". The factor 5/6 mentioned earlier is included so that g'-^^ 
aX Nc = 'i exactly corresponds to the usual nucleon gA-, which has the value 1.2701 ± 0.0025 
[49]. This definition of the axial couplings is convenient in the context of the l/N^ expansion, 
as the differences between the different axial couplings are 0{1/N'^). 



/ \ / \ / \ V 

J — s — I — > — — . — ' 1 — e- — ® — I 1 — . — — 



(1) (2) (3) (4) 

FIG. 2: Diagrams contributing to the 1-loop corrections to the axial-currents. The crossed circle 
denotes the axial-current operator. 

The determination of the axial couplings to require the calculation of the 1-loop 

corrections to the axial current. Only the contributions with no pion pole are necessary, and 
they are given by the diagrams in Fig. IV. The resulting 1-loop contribution to the axial 
currents reads: 

6A\%,^^ = Mr_,,,,(l) + Mr_,oop(2 + 3) + SA\%,,^{A) ., (19) 
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where SA'^^_i^^p{2 + 3) is given by a factor 1/2 times the no-baryon-pole contributions of 
diagrams (2+3). The different contributions read as follows, where one needs to take the 
limits p°,p'° ^ 0: 

Ja 1 



loopK J F3d-1 



^—^ " " — p'^ — Sirin + 5m. 



n,n' 



5A'^_i^^.p{2 + 3) = ^ [C'SZi^ioop + SZi^ioopG"") 



2 

5^r-.oop(4) = -^A(M.)G-. (20) 



Obviously, G"* and SZi^ioop do not commute in general. The pion tadpole integral in the 
last term is given by: 



A(M.) = -^(A,-log-^). (21) 



Notice that the contribution by diagram (4) is actually C(^'^), and thus beyond the degree 
of accuracy of the present calculation. It can serve however as a measure of the size of the 
NNNLO corrections. 

The corrections to the axial currents must scale as 0{NI^) with u < 1. While diagram 
(4) is 0{N^) and therefore consistent in itself, diagrams (1) and (2+3) above are 0{N^). 
As shown in Ref. [28], the offending terms cancel upon adding the diagrams. To test 
the cancellation it is sufficient to take the large Nc limit at fixed M^r. A straightforward 
evaluation leads to: 

SA\%^opmSA\%^^{2+3) 1^ = -^||^ {l[[G^\ G^-^]^/i_z..,(p°, + ■ ■ ■ | , 

(22) 

where ■ ■ ■ indicate further terms which are consistent with the power counting of the 
axial currents. The suppression of the 0{N^) terms is direct consequence of the appearance 
of the commutator of two generators G, which is (9(A^^), when the diagrams are added up. 
In consequence the displayed terms are 0{N^). 

The UV divergent contributions of the individual diagrams read (using the notation 
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-S" = {S{S + 1))''): 



5Af_ 



\UV 



loop\ 



K 9a 
487r2 F3 



+ 



Mr_,„„,(2 + 3) 



uv 



: {^^^ (^QjbQ^a \]^Q3h^g2]^ ^ g2\^ ^ J^2^ [^^G^'']] C'C'^ 

S\ G^''] G'" [G^■^ 5'] ) } 
^^{3M2{G'»,G2} 

G"'G'\&\ S\S^] + [^^ [^^ G^'^jG^^'G'") } 



967r2 F2 



uv 



K 9a 



Mi G' 



(23) 



One notices that only the terms proportional to M| in diagrams (1) and (2+3) diverge as 
proportional to A^^, while the terms proportional to C^p arc 0{N^). Thus only the 0{N^) 
of the terms proportional to M^ need to cancel to give consistency. One can easily check 
that such a cancellation indeed occurs, leaving only terms 0{N^). An explicit evaluation of 
these UV divergent terms using the results from Appendix D finally gives: 



uv 



loop 



'-hf9a 



+ (^4 - 2Ar^(4 + Ar^)G''^ - 7[S^, [S^, G'"]] + 4{^^ G"}) } (24) 

The terms in the Lagrangian needed to renormalize the axial currents are then the fol- 
lowing: 



Bt 



Ui 



These are all the terms which can contribute to the axial currents up to 0(0) which will 
determine the axial couphngs up to 0(^^), i.e., NNLO, which is what is needed for our 
purpose. There are several very important observations concerning the ^-power counting. 
The corrections to the axial couplings start at 0(0) ^"^^ the individual contributions of the 
different baryons in the loop diagrams are also 0[^. Even the difference of different axial 
couplings g^' — 9a"^" starts at 0{^). These latter differences are UV finite. The large 
Nc cancellations do not seem manifest. However, at Nc — 3, where the O^xpansion is used, 
cancellations do occur numerically as shown by Fig. 4 in Section V. Thus, the smallness of 
0{^) terms in the axial couplings is the incipient manifestation of the cancellations in the 
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large Nc limit. If one would consider the strict large Nc limit, the one loop corrections and 
counter-terms considered give the following l/N^ power counting: 

jr-sr" = o(^)+o(^). (26) 

where, as expected, the latter differences are UV finite as in the ^ expansion. 

While in next section a discussion of the nucleon's qa in the context of LQCD results, 
one can make an estimate of the spin-flavor symmetry breaking terms in the axial couplings 
g^^ vs g^'^ using the result for the A width: 

Ta^^n = ^ (^^)' ((m^ - - Mlfl^ (27) 

with the experimental value V i^^T^N{,Exp) = 116 — 120 MeV [49], one obtains g^^ = 1.235 ± 
0.011, which is remarkably close to g^^ = 1.2701 ± 0.0025 [49]. 



V. ANALYSIS OF LATTICE QCD RESULTS FOR BARYON MASSES AND THE 
NUCLEON'S AXIAL COUPLING 

As an application of the present framework of the ^-expansion, this section presents an 
analysis of LQCD results for baryon masses and the nucleon's axial coupling. 

Lattice QCD calculations of the non-strange ground state baryon masses (both of and 
A baryons) have opened the possibility of determining the quark mass dependencies, and 
similarly for the axial coupling of the nucleon. These calculations represent a very fruitful 
ground of applications for ChPT, allowing in particular for a study of the convergence of 
the low energy expansion. Modern dynamical three-flavor {Nf = 2 + 1) calculations of the 
hadron spectrum, and in particular of baryon masses, with fixed strange quark mass and 
variable rriu = md [34-39] are achieving remarkably accurate results in a range of quark 
masses where extrapolations to the physical limit are now possible using effective theory. 
All calculations present similar results for the and A masses, namely, roughly linear 
dependencies of the masses as a function of M^^, and extrapolations to the correct physical 
value within a few percent. For the nucleon axial coupling g^^ the results are particularly 
interesting [40-45] because it shows small dependence in a broad range of M^^. The most 
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recent LQCD calculations for Nf = 2 [40, 42, 45] and Nj = 2 + 1 [41, 43, 44], all agree 
on that observation. An open issue is that all calculations give an underestimation for the 
value of g^^ of about 12% below the experimental value. 

Effects due to finite volume of the lattice have been studied for the observables considered 
here. Those effects are determined primarily by the value of the product LMt^, where L is the 
length of the lattice. For the baryon masses, the rule LM^r ~ 4 [33] seems to be sufficient for 
the volume effects to be negligibly small. On the other hand, for gj^ the LQCD understanding 
of the finite volume effects is not yet complete. According to Ref. [43], gA clearly exhibits 
scaling in LM^ and in lattices with LM^ ~ 4 — 5 the effect on (7^ is a 9 % reduction in 
calculations with 2 + 1 flavors of domain wall fermions and a 25 % reduction in calculations 
with 2 flavors of Wilson fermions. This has led to the current view that LMt^ > 5 — 6 or 
even higher may in fact be needed to reliably determine gA- Finite- volume effects for mases 
and the nucleon axial coupling have been studied in effective theories [50-59]. A detailed 
study of these effects in the present formalism is beyond the scope of this work, and will be 
presented elsewhere [60]. 

In the following, combined flts to LQCD results for and A masses and the nucleon 
gA as functions of Mj,- are carried out. Several studies of the available results have been 
carried out, but only the analysis of the more recent ones are discussed here. For the 
A^ and A masses the results used are those from the PACS-CS collaboration of Ref. [36] 
and the LHP collaboration of Ref. [35]. For gA the results used are those from the LHP 
collaboration [41] and from the ETM collaboration [42]. All collaborations obtain results 
satisfying the constraint LM^ > 4 and for quark masses reaching down close to the physical 
point, in particular for the baryon masses. The flts are carried out only including results 
where LM^r > 4. The analysis of these LQCD results is carried out up to 0{^^) for the 
masses and for gA- The set of Lagrangian counter-terms is the one displayed in 

Eqs. (16) and (25), which are summarized by the following equation: 



5S^^(p° = 0)(5) 




TT 1 



+ 




(28) 
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There are several LECs which to be determined require knowledge of results at different 
values of Nc- With the LQCD results at fixed Nc — 3, those LECs combine with existing ones 
at lower order, making their determination impossible. In addition, because LQCD results 
on g^'^ are not analyzed and the lack of results for g^^, LECs which split the values of the 
different gA^s cannot be fixed either. For instance, the LECs nii and Wi give the sub- leading 
Nc dependence mo, and therefore at fixed Nc — 3 they are absorbed into the fitted value of 
mo- The same will happen with Chfi and W2 with Chf- The LEC /ii is a correction to the 
LO cr-term LEC Ci. Similarly, one can not separate Cq from qa- Therefore, without loss of 
generality, the redundant LECs can be set to vanish. In addition, since the current fits only 
involve the nuclcon's axial coupling, (7^^, not all LECs affecting the axial currents can be 
determined as mentioned earlier. In particular, counter-terms with commutators in Eq. (25) 
only appear in Qa^. Of course, depending on the order in the ^-expansion, the number of 
LECs varies. Specifically, at leading order (LO), that is 0{^) for the mass and 0(^°) for 
the axial couphng, the LECs are mo, Qa, Chf and Ci, at NLO the additional LECs Chfi, 
jii and Cq appear. Finally, at NNLO, that is 0{^^) for the mass and C(^^) for the axial 
coupling, the additional LECs 1x2, zi and Cj^, which are fitted, and wi, W2 and (^^3,4 
that cannot be determined, make their appearance. 

The combined fits to N and A masses and to g^^ up to NNLO for the four possible 
combinations of LQCD results from the collaborations considered here are presented in 
table I, which shows the values for LECs obtained from the fits and the extrapolated values 
for m^v, ?7iA and gA to the physical point. To estimate the theoretical errors, the original 
lattice results are bootstrapped by Montecarlo, and the errors correspond to a 68% confidence 
interval. The LO fits arc restricted to pion masses in the range M^r ~ 140 — 300 MeV, NLO 
fits to ~ 140 - 500 MeV, and NNLO fits to ~ 140 - 600 MeV. 

The following remarks on the fits are in order: 

1. All fitted LECs are of natural size when the renormalization scale is taken to be 
II ~ rup. 

2. Parameters appearing at lower orders, namely mo, gA-, Chf remain stable at higher 

orders, except Ci which changes by more than the estimated 30% when increasing the 
order in ^ of the fit by one unit. 

3. For baryon masses, LQCD data and physical point values are consistent even at LO 
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FIG. 3: Combined fits to PACS-CS [36] and LHP [41] corresponding to the results shown in the 
first row of Table I. LO (black long-dashed line), NLO (black short-dashed line) and NNLO for N 
(blue solid line) and A (red solid line). The bands correspond to the theoretical 68% confidence 
interval. 

where only with three parameters one can extrapolate to the correct experimental 
values. This is not surprising taking into account that a single straight line can be use 
for consistent fit [35]. 

4. For the case of the axial current, cancellations of large contributions from individual 
loop diagrams are very pronounced and the almost flat behavior of qa as a function of 

obtained in LQCD is naturally explained. This is shown in the left panel of Fig. 4 
which shows the finite one-loop contributions to g^^ from each diagram [fi = 700 
MeV). As stated in Eq. (22) this cancellation is exact in the large limit. However, 
at Nc = 3 this cancellation is not exact but still quite pronounced (black curve in 
Fig. 4), and plays the key role in explaining the small dependence in M^. A similar 
cancellation occurs between the contributions of N and A in the loop. This is shown 
in the right panel of Fig. 4. 

5. The physical g^^ cannot be fitted along with the lattice results, instead the lattice 
results and the expansion to NNLO extrapolate to a value 12% smaller value than the 
physical one. The recent LQCD results [61] which reach further down in M^r continue 
that trend. This seems to be therefore an issue with the LQCD calculations rather 
than the effective theory. In fact, it would be very unnatural for the effective theory 

20 



100 



1 


total — 




(1) — 


^^^^ 


(2+3) - — 










1 





200 



300 400 
Mn [MeV] 



500 



600 



1 

0.8 

S 0.6 
'c 

Li_ 

- 0.4 
z 

i 0-2 




-0.2 



N + A in loop ■ 
Only N In loop 



100 



200 



300 400 
Mn [MeV] 



500 



600 




FIG. 4: Finite parts of the one-loop contributions to • The upper left panel shows the individual 
contributions of the diagrams in Fig. 2 up to O{0\ The right panel shows the corresponding effect 
of switching off the contribution of the A in the loops. The third panel shows the effect of removing 
the contributions of the counter-terms to the masses. Throughout [i = 700 MeV. 

to describe a relatively fast changing pattern at low M^r, to continue with the slow 
changing behavior at larger M^r, which it can describe in a natural fashion. 

6. A fit restricted only to masses gives too small a value for (7^^, namely, qa ~ 0.5 — 0.8. 
A realistic value can only be obtained with the combined fit. 

7. Predictions for ^f^^ and g^^ cannot be made without the corresponding LQCD results. 
However, the results at the physical point from Eq. (27) suggest that these are going 
to be very similar in value to g^^- Further efforts to study this coupUngs in LQCD 
will be very useful. 
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8. In the masses one finds that above > 350 — 400 MeV there is a significant can- 
cellation between the contributions of the one-loop diagram and the counter-terms as 
shown in Fig 4, which must be taken as an indicator of the range of convergence of 
the expansion. The masses counter-terms are 0{M^N^) = 0{C,^). 

9. Using the results of the fit at one obtains for the a terms evaluated at the 
physical pion mass: 

= 65 ± 8MeV 

= 90± 9MeV. (29) 

10. It must be emphasized that the results obtained here are very similar to those ob- 
tained in works where the A has been included explicitly [16, 19, 20, 62-64]. The 
main advantage of the present approach of the ^-expansion is its systematic character, 
which in particular will be more prominently shown when carrying out higher order 
calculations than the ones considered here. 

VI. DISCUSSION AND CONCLUSIONS 

Chiral symmetry and the large limit are of fundamental conceptual importance in 
QCD. The former is known to play a crucial role in light hadrons, and there are multi- 
ple indications that the latter is also important, in particular for baryons. It is therefore 
very important to have a theoretical framework where both of these aspects of QCD are 
consistently incorporated. This is possible with the combined 1/Nc and Chiral expansions 
of QCD, which in the baryon sector is implemented with the effective theory discussed in 
this work. A particular power counting, the ^-expansion, which links the l/N^ and low 
energy expansions as 1/Nc = 0{C,) = 0{p) is proposed as the most realistic one for studying 
baryons at Nc = 3. Results for the masses and axial couplings at NNLO have been given, 
and applied to current LQCD results. 

The (^-expansion at NNLO clearly provides a satisfactory description of the LQCD results, 
and in particular it illuminates the mild dependence of the axial couplings on the quark 
masses as a result of important cancellations, which had been realized in various previous 
analysis by various groups. It is important to complete the study in SU{3), in particular 
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because the one-loop contributions to the baryon masses become larger in magnitude, and 
a smaller range of convergence is expected. These results will be presented elsewhere [65] 

The deficit in at the physical point seems to be a LQCD issue rather than a problem 
of convergence of the effective theory. In particular because the ^-expansion is especially 
well behaved for qa- It is however a puzzling problem, because the likely candidate for the 
explanation, namely finite volume effects, is expected to be significantly somewhat than the 
effect observed. 

In addition to the tests LQCD can provide on the quark mass dependencies, it is the 
ideal tool to test the N^. behavior of QCD. Baryon LQCD is becoming accessible at varying 
values of [66], which is a promising development. 
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Appendix A: Spin-flavor Algebra 



The 4Nj — 1 generators of the spin- flavor group SU{2Nf) consist of the three spin gener- 
ators S**, the Nj — 1 flavor SU{Nf) generators T°-, and the remaining 3(A^| — 1) spin/flavor 
generators G*". The commutation relations are: 

[S\ S^] = ie,,uS\ [^^ T'] = tU.T'^, S'] = 
[S\ G^"] = ze.jkG^'^, [T«, G''] = ifabcG^' 

^Qta^ Qjb^^ ^ t_^ijjabcrpc ^ _l_^ab^tjkgk ^ t_^ijk ^abcQkc (^^^^ 

For two flavors one has the isospin generators a = 1, 2, 3. 

In representations with indices (baryons), the generators G**^ have matrix elements 
0{Nc) on states with S = O^N^!). A contracted 5'f/(4) algebra is defined by the generators 
{S^,I°',X'^°'}, where X*" = G^°'/Nc. In large Nc, the generators X*" become semiclassical as 
[X^'^jX-'^] = 0{1/N^), while having matrix elements 0{1) in baryon representations. 

Appendix B: Non-linear realization of chiral symmetry and spin-flavor transforma- 
tions 

In the symmetric representations of SU{4) the baryon spin-flavor multiplet consists of the 
baryon states with I = S. In particular, isospin transformation will act on the spin-flavor 
multiplet in an obvious way on those states. This permits a straighforward implementation of 
the non-linear realization of chiral SUlI^) x SUr{2) on the spin-flavor multiplet. Defining as 
usual the Goldstone Boson fields tt'^ through the unitary parametrization u = exp(i^^^), for 
any isospin representation one defines a non-linear realization of chiral symmetry according 
to [3, 4]: 

(L, R) ■.u = u = Ruh\L, R, u) = h{L, R, u)uL\ (B1) 

where {L,R) is any SUl{2) x SUii(2) transformation. This equation defines h, and since h 
is an isospin SU{2) transformation itself, it can be written a.s h = exp(zc"J"). The chiral 
transformation on the baryon multiplet B is then given by: 

(L, R) :B = B' = h{L, R, u)B. (B2) 

On the other hand, spin-flavor transformations of interest are the contracted ones, namely 
those generated by {^\J'^,X^" = ^^G'"}. While the isospin transformations act on the 
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pion fields in the usual way, and the spin transformations must be performed along with 
the corresponding spatial rotations, and the transformations generated by X" are defined 
to only act on the baryons. 

Appendix C: Tools for building effective Lagrangians 

The effective baryon Lagrangian can be expressed in the usual way as a series of terms 
which are SUl{'2) x SUr{2) invariant (upon introduction of appropriate sources; see for 
instance [67] for details). In addition, implemented in the effective Lagrangian is the ap- 
proximate SU (4) symmetry and its breaking as a power series in 1/Nc [27]. The fields in the 
effective Lagrangian are the Goldstone Bosons parametrized by the unitary SU{2) matrix 
field u and the baryons given by the symmetric SU{A) multiplet B of / = S* fields. 

The building blocks for the effective theory consist of low energy operators, and spin-flavor 
operators. 

The low energy operators are the usual ones, namely: 

X = 2Bo{s + ip), x± = u'^xu'' ±ux^u 
p^y = df't - d^t - z[r, r], Fj^'' = d^^r" - d^r^" - z[r^, r"] (CI) 

where is the chiral covariant derivative, s and p are scalar and pseudo-scalar sources, x± = 
2M^ + ■ ■ ■ , and £^ and are gauge sources. The spin-flavor operators are tensor operators 
consisting of products of the spin-flavor generators. These operators can be reduced by 
means of the commutation relations to forms which only contain anti-commutators. A set 
of identities shown in Table nnnn permits one to arrive at sets of basis operators at each order 
in 1/Nc for a given spin/isospin tensor type of operator. The l/N^. order v of an operator 
reduced as mentioned is z/ = — 1 — k, where n is the number of generators appearing as 
factors in the operator (one then says that the operator is an n-body operator), and k is the 
number of generators G*" in the product. 

The leading order equations of motion can be used in the construction of the higher order 
terms, namely, iD^B = {^S{S + 1) + f A'cX+)B, and V^u'^ = |x- 
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Appendix D: Matrix elements of spin-flavor operators in the symmetric represen- 
tations of SU(4) 



The evaluation of the matrix elements of spin-flavor operators in the present work can be 
carried out starting from the following matrix elements of the spin-flavor generators in the 
totally symmetric representation of SU (4) corresponding to the Young tableux with a single 
row of Nc boxes. The basis states of the symmetric representation consists of the states with 
I = S, namely | 5* 6*3/3), where S3 and I3 are the spin and isospin projections respectively. 



{S'S'31'3 I G^'^ I 553/3) = -y^^7^C(iVc,5,5')(5 53,H I 5'5^)(5/3,la | 57^), (Dl) 



where C(A^c, S, S') = y/{2 + N,)^ - (S - S')'^{S + S' + 1)2 [68]. The products of generators 



can be reduced by means of the use of the commutation relations, and further, for matrix 
elements in the symmetric representation, via the reduction rules [69], which for convenience 
are displayed in the following Table II. 

TABLE II: SU (4) operator identities in the totally symmetric irreducible representation {Nc, 0, 0) 
of 5^7(4). The last column gives the operator's quantum numbers {J, I) under SU{2) x SU(2) 



I s' I 553/3) 



VS{S + T)5ss'Sj,i;^{SS3,li I S'S',) 



{S'S'^l'^lPlSSsh) 



y/S{S + l)5ss'Ss,s;,{Sh,la \ S'l'^) 




{S\S'} - {/",/"} = 



(0,0) 



{5% S'} + {I", I"} + 4{G^'^, G'"} 



|iVe(4 + iVc) 



(0,0) 



2{S\G'''} = (2 + iVc)/' 



a 



(0,1) 



2{ja^(^ia| ^ (^2 + Nc)S' 

^{5^ r} - e^i'=e«f'c{G'^«, = (2 + Nc)G' 
4{G-,G^^}|/=2 = {/^/''}|/=2 



(1,0) 



(1,1) 



(1,1) 



(0,2) 



4{G-,G^■nlJ=2 = {S^5n|J=2 



(2,0) 
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1. Useful matrix elements 



It is always convenient to express matrix elements in terms of reduced matrix elements 
(RMEs) defined in the ordinary Wigner-Eckart fashion [70]. The RMEs defined here are 
with respect to SUspm{'2) x SUi{2). For matrix elements in the symmetric representation of 
spin flavor the Wigner-Eckart theorem reads: 

{S'S!,I', I I ^^3/3) = ^^^2g^|i^^ ^^3,JJ3 I ^'^^)(^/3,Il3 I S'l',), (D2) 

where O is an SUs{2) x SUj{2) irreducible tensor operator, and (S" || Oj || S) is the 
reduced matrix element. Note that the notation || S) indicates the spin-flavor states in the 
symmetric representation {Nc, 0, 0) with I = S. The reduced matrix elements of the SU (4) 
generators read: 



{S' 1 


1 s 1 


\S) 


{S' 


1 1 1 


\S) 


{S' 1 


1 G 1 


\S) 



6{s,s',i}^Vi^S+l){2S' + l)aNc,S,S'), (D3) 

where 5{5,s",i} = 1 if | S* — S" |< 1 and otherwise vanishes. 

Reduced matrix elements of the operators involving the projects Vn are easily obtained 
using that P„ = Z]s„3,/„3 I Sn, Sn3ln3){Sn, Sns^ns \, and the SU{2) re-coupling resuhs [70]. 
For the masses the relevant such RME becomes: 

{S II G^'^VnG"' \\S) = ^^{Sn II G II 5)2. (D4) 

For the axial currents the following RME is needed, namely: 



{S' II G^^Vn'C'VnG^'' II S) 



_ 2 

S Sfi 1 

Sn' S' 1 

X {S' II G II Sn'){Sn' II G II Sn){Sn \\ G \\ S) . (D5) 



Various reduced matrix elements which appear in the evaluation of the UV divergent 
pieces of the one-loop contributions to the self energy and to the axial currents are given 
below. They are obtained using the results given in the above Eqs. (D4) and (D5): 

{S'\\G'\\S) = Uss'i2S + l)i-S{S + l) + ^iV,(4 + iV,)), (D6) 
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where = G'^G'". 

{S' 1 1 GS^G \\S) = ^dssi2S+l) (^N,{A + N,) + ^(-40 + + 3N'^)S{S + 1) - {8(8 + 1))^ ) , 



(D7) 



+ (-7+lN^+^N'^ms+i)r-l{s(s+i)ry (ds) 

In the foUowing, 8' = 8 or 8 ± 1. With obvious notation: 

TV , 

{S' II \\S)^ 5ss'{l + -j){2S + lWSiS + l), (D9) 

obtained using the corresponding reduction relation in Table II. 

{S' II S'P II S)^6ss'{'2S + l)S{S + l), (DIO) 

{8' II II 5) = -{8{8 + 1) + + 1) - 2)(5' II G II 5), (DU) 

/^/ II II 8) = ^{8' II G II 5) 

16 

X (3iVe(4 + Ar,) -4(2 + 5(5+1) + + (D12) 

1 /S 
(5' II GG"'8^G II 5) = -(5' II G II 5) i-Nc{4: + Nc) 

+ 5(5 + l)(-5(5 + 1) - 5 + In,{4 + iVj + J5'(5' + 1)) 

o z 

- S'{S' + 1)(1 + ^5'(5' + 1)^ , (D13) 

(5' II GG'^S^G II 5) = ^(5' II G \\ S) 

X Q7V,(4 + TV,) + (5(5 + 1))'(-16 + ^7V,(4 + TV,) - 35'(5' + 1)) 

+ (5'(5' + l))2(l-^5'(5' + l)) 

+ 5(5 + 1)(-12 + ^7V,(4 + AT,) + Z(5'(5' + 1))^)^ , (D14) 

(5' II GS'^G'^'S'^G II 5) = ^(5' II G II 5) 

X (^(2 + Ar,)2 + ^(-16 + 57V,(4 + iVc))(^('5 + 1) + 5'(5' + 1)) 

- (5(5 + 1))2(9 + 5(5 + 1)) - (5'(5' + 1))2(9 + 8\8' + 1)) 

+ ^(-16 + Ar,(4 + Ar,))5(5 + l)5'(5' + l)y (D15) 
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Finally, using that for any spin and isospin singlet operator (not necessarily an SU (4) singlet) 
Q, {S' II QO^ \ \S)^ 2^i^' II 2 II '^')('^' II Oi II one can easily obtain the rest of the 
matrix elements involved in the calculation of the axial currents. 
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